Homodyne detection of matter-wave fields 
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A scheme is proposed, that allows one for performing homodyne detection of the matter-wave field 
of ultracold bosonic atoms. It is based on a pump-probe lasers setup, that both illuminates a Bose- 
Einstein condensate, acting as reference system, and a second ultracold gas, composed by the same 
atoms but in a quantum phase to determine. Photon scattering outcouples atoms from both systems, 
which then propagate freely. Under appropriate conditions, when the same photon can either be 
scattered by the Bose-Einstein condensate or by the other quantum gas, both flux of outcoupled 
atoms and scattered photons exhibit oscillations, whose amplitude is proportional to the condensate 
fraction of the quantum gas. This setup allows one, for instance, to perform thermometry of a 
condensate or to monitor the Mott-insulator/superfluid phase transition in optical lattices, and can 
be extended in order to measure the first-order correlation function of a quantum gas. The dynamics 
here discussed make use of the entanglement between atoms and photons, which is established by 
the scattering process, in order to access detailed information on the quantum state of matter. 



I. INTRODUCTION 
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FIG. 1: Setup for performing homodyne detection of matter- 
wave fields. A gas of identical, bosonic atoms is confined 
in two distinguishable regions of space. In the left region it 
forms a Bose-Einstein condensate, that acts as reference sys- 
tem for measuring the condensate fraction of the right system. 
A pump-probe laser setup outcouples the atoms from both 
wells. The photon flux of the probe beam is constituted by 
the photon scattered in the outcoupling process and exhibits 
oscillations whose amplitude depend on the condensate frac- 
tion of the right system, and which vanish when this is zero. 
The first-order correlation function of a quantum gas can be 
determined in an extension of this setup, where the shaded 
regions denoted hy xpL and tpR correspond to the spatial re- 
gions of the quantum gas trapped in a single well. In this 
case, the amplitude of the photon flux oscillations provides a 
measurement of the spatial correlation function as a function 
of the distance d. Detail on the parameters are reported in 
the text. 

The measurement technique, which is mostly employed 
in experiments with ultracold atoms, is based on time of 
flight [H, Used in different setups, it allows one to de- 



termine the quantum state of matter . The most re- 
cent demonstration of in-situ imaging of atoms in optical 
lattices 0] has opened further exciting directions for 
probing and manipulating atoms in periodic arrays Q . 
In addition, the remarkable progress in coupling ultracold 
atomic gases with high-finesse optical resonators [lol - [l3j 
has generated a renewed interest in revealing the corre- 
lation functions of matter by photodetection of the scat- 
tered light. This progress opens, amongst several others, 
the possibility to monitor properties of the quantum state 
of matter in a non-destructive way p^| - |T7| . 



Photons emitted by Rayleigh scattering, and in general 
in a pump-probe type of experiment, have been used in 
order to characterize the quantum state of atomic gases 
p"8l - [2^ . They deliver information on the structure form 
factor of the quantum gas and thus on the density and 
density-density correlations 0, [22| - [23 | . Proposals for op- 
tical detection of certain quantum states of matter have 
appeared in the literature 23 27j. To our knowledge, 
however, a setup has not yet been identified, which would 
allow one measure the mean value of the atomic field op- 
erator by means of the scattered light. This measure 
would be analogous to homodyne detection of light fields 
[281] . and specifically allow one to determine the conden- 
sate fraction of a quantum gas [2^ by photodetection. 

In this article we propose a setup in which detection of 
the photons scattered by a quantum gas allows one to de- 
termine the mean value of the atomic field. The setup for 
determining the condensate fraction is sketched in Fig. [TJ 
and is based on an interferometer for Bose-Einstein Con- 
densates (BEC) realized in Refs. [sO, HH. We establish a 
direct connection between our theoretical model and the 
experimental setup in Ref. j30l , IsTj . and then focus on 
a scheme in which the second system is not necessarily 
a BEC. We then show that this setup allows one to de- 
termine the condensate fraction of the second system by 
measuring the flux of the scattered photon. Furthermore, 
we argue that this scheme can be extended to determine 
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the first-order correlation function of an atomic gas, when 
the two illuminated regions belong to the same quantum 
gas. The dynamics is based on first creating entangle- 
ment between the scattered photon and the scattering 
atom of the quantum gas, and then on performing an op- 
eration similar to a quantum eraser [s^ , thereby allowing 
one to measure the correlation functions of matter in the 
flux of the scattered light. 

This article is organized as follows. In Sec. HI] the setup 
and the theoretical model are introduced. In Sec. IIIII a 
formal connection is derived between the photon flux and 
the correlation functions of the scattering atoms. In Sec. 
llVl the photon flux is evaluated for three exemplary cases: 
when the scattering systems are two Bose-Einstein con- 
densates at zero temperature, when they are at different 
temperatures, and when one of the systems is composed 
by ultracold atoms in a two-dimensional optical lattice. 
In Sec. |V] it is discussed how this setup may allow one 
to measure the first-order correlation function. The con- 
clusions are drawn in Sec. I VII The appendices report 
details of the calculations presented in Sees. IIIII IIVI and 
lYl 



II. THE MODEL 

A realization of the setup we consider in this paper is 
depicted in Fig.[T] Here, ultracold bosonic atoms of mass 
m are initially prepared in the stable electronic state |1) 
and confined by the (state-dependent) potential Vi(r). 
The atoms, that are illuminated by the lasers, are coher- 
ently coupled to the electronic state |2), which is stable 
and not confined: the outcoupled atoms propagate freely 
in space. 

In this section we introduce the Hamiltonian and the 
physical quantities which are at the basis of our analysis. 



A. Hamiltonian 

We denote by H the Hamiltonian governing the dy- 
namics of photons and atoms, which we decompose into 
the sum 



cmf 



int 



shift : 



(1) 



with Hq the Hamiltonian governing the atoms dynamics 
in absence of interaction with the lasers, Hcmf the Hamil- 
tonian for the free, transverse electromagnetic (e.m.) 
field, while Hint describes the coupling between states |1) 
and 1 2) induced by coherent Raman scattering. Hamilto- 
nian -ffshift includes the dynamical Stark shift due to off- 
resonant coupling between atoms and e.m. -field modes. 
These Hamiltonians are reported in a second-quantized 
description of matter- and photon-fields. In detail. 



gives the energy of the transverse electromagnetic field in 
free space (and is reported without the vacuum energy) , 
where A labels a mode of the e.m.-field at wave vector k, 
polarization e _L k, and frequency Wk — c|k|, with c the 
velocity of light. The operators a a and a\ annihilate and 
create, respectively, a photon in mode A, and obey the 
bosonic commutation relation [aA,a]^,] = i5a,a'- 

The atomic Hamiltonian Hq is conveniently rewritten 
as Ho = Hi + H2 + Hi2, where 



H 



0=1,2} 



dr-^Wr) 



2m 



+ 1 j drV.](r)V'](r)V',(r)^,(r) (2) 

describes the atoms dynamics when the atoms are in 
the electronic state \j) (at frequency ujaj)- Here V'j(r) 
and V'j(j') the annihilation and creation operator, 
respectively for a bosonic atom at position r and in 
the electronic state \j — 1,2), obeying the commuta- 
tion relation [tpj{r),-ipl{r')] — Sj^kS{r — r'). Parameter 
gj — ATth^asj /m is the strength of s-wave scattering be- 
tween the atoms in state with the s-wave scat- 
tering length in the corresponding electronic state poj . 
Hamiltonian term H12 describes s-wave scattering be- 
tween an atom in electronic state |1) and an atom in 
electronic state |2): We do not report its explicit form as 
we will consider situations for which these types of col- 
lisions can be neglected. The interested reader can find 
the corresponding term, for instance, in Refs. ,33.. 34] . 

The Hamiltonian terms describing the interactions be- 
tween photons and atoms are given using normal order- 
ing [35] . They include a term describing coherent Raman 
coupling between the states ]1) and ]2), in which a pho- 
ton in mode 1 is scattered into mode 2 and vice versa. 
Raman transitions follow from a pump-probe type of ex- 
citation, and the corresponding Hamiltonian reads 



Hint dr 7(r)4'02(i")V^i(i")ai + H.c 



(3) 



with the position-dependent coupling strength 7(r). This 
interaction can be tailored by means of a laser, pumping 
the atoms, and a cavity mode, acting as a probe, that are 
both far-detuned from the atomic excited state but are 
set close to resonance with an atomic Raman transition. 
For later convenience we denote by the quantity 

UJ12 — LOl ~ LO2 

the difference between the frequencies of the two e.m.f.- 
modes. Moreover, we assume a;i2 ^ uia2 ~ Wai, i.e., the 
two lasers drive quasi-resonantly the Raman transition 
coupling the stable electronic states ]1) and ]2). The 
interaction with the pump and probe laser also induces a 
dynamical Stark shift, whose corresponding Hamiltonian 
reads (in normally-ordered form) 



H 



shift 



dr7j(r)a]V'](r)V'j(r)aj , (4) 
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where the parameter 7^ (r) has the dimension of an angu- 
lar frequency. It can be written as 7j(r) = 2|rij(r)p/A, 
with rij(r) the Rabi frequency of the mode couphng 
the transition |j) — ?> |e) and A the detuning between 
the mode and the atomic transition frequencies (the 
Lamb shifts are included in the frequencies of the atomic 
states). 

In the following we will assume that the pump and 
probe fields are traveling waves with wave vectors kj . In 
this case, the Raman- coupling strength 7(r) reads 

7(r) = 7oe''i-"- (5) 
with q = ki - ka and 70 = 217i(r)n2(r)*/^- 



B. Spin-dependent potential 

We will now provide more details on the potential 
which confines the atoms. We will assume that the atoms 
in state 1 2) propagate freely, namely, V2 (r) = const in Eq. 
([2]), which we set equal to zero. 

The atoms in state |1) are trapped by potential Vi(r). 
This potential confines the atoms in two spatially sepa- 
rated regions around the two potential minima, which are 
at distance d and, specifically, are localized at the points 
tl = —{d/2)x and tr — {d/2)x. The potential can be 
then decomposed in the sum Vi (r) Vl (r) + Vr (r) , with 
Vj{r) the potential center at rj and j = L,R. The two 
atomic clouds at each well are initially uncorrelated, and 
there is no tunneling between the two spatial region. 

In principle, hence, the atoms in the two regions of 
space are distinguishable. It is then useful to consider 
the partition -01(1") = ?AL(r) where we denote by 

V'L(r,i) = i}i{r, t)0{-z) and ■ipR{r,t) = Vi(r,i)6l(z) the 
field operators which do not vanish on the left and right 



region of space, respectively, such that [■i/'j(r), ■0^(1"')] = 
dj^kS{r — r') for r, r' ^ [36|. Using these definitions, we 
can write Hi = Hl + Hr, with 



+ y,(r) 0,(r) 



+ f / drV4(r)^Ur)^,(r)V'j(r), (6) 



and where we have set LUai — 0. With this representation, 
the Hamiltonian describing the interaction with the lasers 
takes the form 



dr 



(7) 

We note that our model is an extension to the one studied 
in Refs. [13, HI], where the authors studied the outcou- 
pling of atoms from a single BEC by means of classical 
Raman lasers. By considering the quantum dynamics of 
the interaction between photons and atoms, we take into 
account the quantum fluctuations of the light field due 
to the scattering process, which are discarded in (37l.[38|. 
We will indeed show that these fluctuations give access 
to some correlation functions of the scattering system. 
C. The scattered field 



We now study the properties of the scattered photons 
by considering the photon field operators in Heisenberg 
picture. The Heisenberg equations of motion are deter- 
mined assuming that the coupling between matter and 
photons is sufficiently weak to be treated in second-order 
perturbation theory. In this limit the operators for the 
photonic modes read 



J 



ai( 



(t) = ai(0)e-'("i*+'^i(*» -ie-'"i* dre''^^^-' / dr 7(r) tl;^°"{r,T)+tl;^°"{r,T) (r, ^)a2(0) , (8) 



a2( 



(t) = a2(0)e-'("^*+*^(*» -ie-'"^* / dre"''^"" / dr 7(r)4°^^(r, r) (r, r) + V^"^ (r, t)1 ai(0) , (9) 



where (r, t) = eyip{iHoT/h)^pj{r,0)eyip{—iHoT/h). 
The physical origin of the individual terms on the Right- 
Hand Side (RHS) can be simply identified. The first 
term on the RHS of both equations is the free-field com- 
ponent. It is characterised by a time-dependent phase 
4>j{t), which is proportional to the atomic density in the 
electronic state When the Born approximation can 
be performed, it is a density-dependent phase shift of the 
field mode, with the form 



t 

<^iW = E /dr /dr 71^^-^ ^(r,T)4°)(r,T) , 

t 

Ut) = JdrJ dr72(r)0f^(r,T)vf (r,T). 
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Since all atoms are initially prepared in the electronic 
state |1), this density-dependent shift can be measured 
in the frequency shift the field-mode 1 after it has in- 
teracted with the atomic medium. This would allow 
a quantum-non-demolition measurement of the density 
of the medium integrated over the path along which 
light propagates Q. It is interesting to note that in 
the considered setup this shift contains also the operator 
ip^^'^\r,T)^^^\r,T) and its adjoint. Hence, its measure- 
ment would reveal tunneling events between the wells. 
This is however not relevant for the case considered in 
this work, since we assume that there is no tunneling 
between the two separate wells. 

The second term on the RHS establishes a direct 
proportionality relation between the photonic and the 
matter- wave fields 02 and V'i(r) (ai and ■(/'2(r))- It shows, 
in particular, that the source term of field ai is the co- 
herent overlap of the field scattered from the right and 
from the left wells. In this shape the outcoupling pro- 
cess in this system resembles a beam-splitter operation. 
This property establishes an analogy with an interfero- 
mctric setup, which we will exploit in order to perform 
homodyne detection of matter-wave field. Differing from 
a simple interferometer, however, atoms and photons are 
correlated by the scattering process. This important dif- 
ference is at the basis of the dynamics we observe. 

Before we discuss the signal at the photodetector, we 
will introduce few approximations which will notably 
simplify the treatment. In first place we neglect atomic 
collisions between outcoupled atoms in state |2), since 
we assume that the gas of outcoupled atoms is at very 
low densities. This regime also permits us to neglecting 
collisions between outcoupled atoms and trapped atoms 
in state |1) [11]. 



III. PHOTON FLUX 

The quantity we analyse in this paper is the flux of 
photons of mode 2, namely, the rate of change of the 
photon number in the mode at frequency U!2- Formally, 
the photon flux is given by the equation 



dt 



al{t)a2{t)) 



(10) 



where the expectation value (•) is taken over the initial 
state of the atoms and the e.m. -field. The photon flux, 
integrated over the detection time, gives the integrated 
intensity of the field at the detector. It can be verified 
that 



^ / dr(^|(r,t)^2(r,i)) 



(11) 



where the proportionality factor is the mean number of 
photons in mode 1. This equality shows indeed that fiux 
of scattered photons and of the corresponding outcoupled 
atoms carry the same information. 



For the specific setup we consider the photon flux can 
be rewritten as the sum of two contributions. 



F = Fb+Fi. 



(12) 



with 



Fsit) = FRe ^ J dr J dr' 

j=R,L 

t 

X J di7(r, t; r',t')G,, (r, t; r', t') , (13a) 


t 

Fi{t) = FRe j dr j dr' j dt' f{r,t;r' ,t') 



X (GLK(r, t; r', t') + GflL(r, t- r', t')) ,(13b) 

where F — 2(a|ai)7Q is a scaling factor, proportional to 
the number of photons in mode 1. Here, 

G,fe(r,t;r',t') = (^](r', OV'fe(r, i)) , (14) 
/(r,t;r',t') = ^Mr-r')-u.Mt-t')] U^(r' ,t')4{r,t)\ 



(15) 

where the initial state is the vacuum state of the e.m. 
field except for mode 1 and 2, which are assumed to be 
in coherent states with non-vanishing photon number, 
while all atoms are in internal state |1) and at equilibrium 
in the grand-canonical ensemble at temperature T. The 
atoms in state |2) propagate freely and do not undergo 
collisions, since we assume that the density is very low. 
Therefore, Eq. ([T5|) can be cast in the form 



/(r,i;r',t')= / ^^^e' 



i[(q-k)-(r-r')-(wi2-Wk)(t-*')] 



with 



(16) 



2 m 



the recoil frequency and 



W12 — ^12 ^ ^a2 

the two-photon detuning. 

In order to obtain explicit expressions for the corre- 
lation functions Gjk{r,t;r',t') it is convenient to work 
in the interaction picture with respect to the grand 
canonical ensemble Kq = Hg — J2j=L r Mi-^j i where 
fij is the chemical potential of the atoms in either the 
right or left cloud and J^j = J dr tpj{r)tpj(r). The 
new atomic field operators are obtained from the ones 
in Heisenberg picture with respect to Hq by replacing 
^j{r,t) — > '0j(r, i)e^'s''^*. From now on we will denote 
by ipj (r, t) the atomic field operators in Heisenberg pic- 
ture with respect to Kq. 
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We now discuss how the considered setup can be used 
in order to measure the mean value of the atomic field 
operator by means of photons. We first note that the pho- 
ton flux is the sum of two terms: the component Fb, that 
is the sum of the flux from each well, and the component 
Fi, that arises from the coherent superposition of a pho- 
ton (atom) scattered by the wells. We will denote Fb by 
"background contribution" and Fj by "interference con- 
tribution". In absence of initial correlations, this latter 
term is proportional to the product of the mean value of 
the field operators {ipL{r,t)){'ipfi{r,t))* . Let us now con- 
sider the case in which, say, the left well confines weakly- 
interacting atoms forming a Bose-Einstein condensate. 
Under this assumption the atomic field operator can be 
written as [l^ 

V'L(r,t) =e-*^^*($L(r)-f <5VL(i-,t)) , (17) 

where 'i?L(r) is the macroscopic wave function, which 
solves the Gross-Pitaevskii equation for the quantum gas 
in the left well of the potential (with chemical potential 
IIl). Then, when the order parameter of the left conden- 
sate is known, the interference term will deliver the mean 
value of the field operator in the right well. 

In the following we use Eq. (|17p in the equations for the 
background contribution to the photon fiux, Eq. (II3ap . 
and for the interference term, Eq. (|13bl) . In particular, 
we will take 



(18) 



and assume that both fii'^) and lpl are real valued. 
Hence, fL{'^)'^ is the density of condensed atoms, while 
the phase i^sl is assumed to be constant in space (hence 
discarding the possibility of superfiuid currents in the 
Bose-Einstein condensate [29|). 



A. Background contribution 

The background contribution can be decomposed into 
the sum of the photon fiux from the left and from the 
right well, Fb = Fl + F^. Since we make a specific 
assumption on the quantum state of the gas trapped in 
the left well, we can provide an explicit form for F^. For 
this purpose, we consider the integrand Gll{y, t; r', t') in 
Fl and observe that it can be splitted into two terms: 



GLL{v,t;r\t')^e 



- Ji-iL(t-t')/h 



accounts for the contribution of the noncondensed atoms 
to the photon fiux. Here, 

oo 

^^,^t(r,r',c.)= / dTe'"-([(5V(r,T),<5V^(r',0)]) 



(19) 



where 



5GLL{v,t;v\t') = {5^l;{{vMLv\t')) 



dw 
2^' 



-\LO(t-t') 



^o(a^)A,^,,^t(r,r',-c^) (20) 



is the spectral density [40|, with Nq{ijj) = [e 



fihuj 



(21) 

1]-^ 



Similar expressions can be found for the background con- 
tribution of the quantum gas in the right well. 

In order to gain insight on the equations we just de- 
rived, we assume that the gas can be considered homo- 
geneous. This assumption allows us to write the spectral 
density as A^^^.^^t (r, r', cj) = A^^^.^^t (|r - r'|, cj). We 

then inspect the equations for the fiux using the momen- 
tum representation for the atomic field operator. The 
terms giving the contribution to the ^^hoton flux from 
the left cloud can be cast in the form 



Fdt) 
SFL{t) 

where A 



27rrRe 
2nTVRe 



dk 

(2^ 
duj 
2^ 



\fL{k)\'d\n-uJw+^), (22a 
dk 



(27r)3 

x^o(^)i5^,5^t(k, -uj)~S\uj 



(22b) 



(k, uj) is the spectral weight function. 



which is the Fourier transform of A 



(r,w) and 

gives strength of the collective excitations at wavevec- 
tor k and frequency w, weighted by the Bose func- 
tion Nq{(jj) [40] • The other parameters are the vol- 
ume V in which the atoms of the left well are confined, 
the Fourier transform of the macroscopic wave function 
/^(k) — J dre~''' ''/L(r), and the frequency 



il = U}i2 — 

n 



(23) 



giving the effective detuning of the laser from the collec- 
tive transition. Moreover, in Eqs. (|22l) we have intro- 
duced the quantity 



(24) 



that is proportional to the diffraction function (5*(a;) — 
sm{ujt/2)/{7TUi), enforcing energy conservation for long 
interaction times t. In particular, S*(u}) — >■ d{u}) for 
t — >■ oo, with S{uj) Dirac-delta function [42|. Definition 
contains a time-dependent phase and highlights that 
the corresponding factor tends to unity in the limit in 
which energy conservation applies. Nevertheless, since 
we are going to also consider intermediate times, we 
will systematically keep this time-dependent phase in the 
equations for the photon flux. 

The integrals in Eqs. ([2^ runs over all values of the 
atomic momentum k. The integrands are the product 
of the momentum distribution at k' — q — k and of 
the diffraction function: The flrst accounts for the 
effect of the photon recoil hq on the atomic distribution 
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due to the outcouphng process, while the diffraction 
function imposes energy conservation in the scattering 
process. Equations ([22k ) and (|22b) show thus that the 
contributions to the flux come from the total number 
of condensed and noncondensed atoms, respectively, 
which fulfill energy and momentum conservation of the 
scattering process. 

We remark that according to Eq. (|TT|) , the flux evalu- 
ated in Eqs. (1^^ gives also the corresponding component 
of the atomic flux. The latter agrees with the correspond- 
ing expressions derived in Ref. 



371 for the atomic flux 



outcoupled from a single BEC by classical flelds. 



Interference contribution 



while GRL{r,t;r',t') = [GLR{r' ,t';r,t)Y 



From Eq. (|25p one observes that these correlation func- 
tions are proportional to the mean value of atomic fleld 
operator (V'_R(r, t)). When this is zero, the interference 
contribution vanishes. Otherwise, the amplitude of this 
term is proportional to the condensate fraction of the 
quantum gas in the right well. We remark that the non- 
condensed atoms in the left well do not contribute to the 
signal, since (i) there are no initial correlations between 
the atoms in the left and right wells, and (ii) the mean 
value ((5V'L(r,t)) = 0. 



In order to determine the interference contribution to 
the photon flux, Fi(t) in Eq. (jl3b[) . one needs the ex- 
plicit form of the correlation functions GLR{^,t;r' ,t'), 
GRL{i^^t\T' ,t'). For the considered setup, however, one 
can already make general statements. In absence of ini- 
tial correlations, in fact, they are the product of the mean 
value of the field operators in each well, and thus take the 
form 

Gifl(r,i;r',t') = e'(^-*-^«*')/^4(r,0)(^fl.(r',t')> 

= e'(^-*-^«*')/V'^-/L(r)(^^(r',t')> 

(25) 



In the shape of Eq. ((25|) the analogy with homodyne 
detection, as it is performed with light fields, can be 
drawn: the condensate in the left well plays the role of 
the local oscillator 28]. Identifying the analon of the 
phase of the local oscillator is however a more delicate 
issue, that deserves some more analysis. For this pur- 
pose, we first assume that the mean value {ipR(r,t)) can 
be evaluated within a mean-field approach, such that 
{ipR{r,t)) = fR{Y)e^f'^, with /ii(r) real-valued function 
and ^pr real constant. The integral in Eq. (|13bl) can then 
be cast in the form 



Fi{t) = 27rrRe 



dk 

(2^ 



/L(k)*/j^(k)^*(f^-C^k+q) 



-i{SfJ.t — ipLR) 



/fl(k)*/L(k)^*(f^-5M-^k+q) ,(26) 



with 



Discussion 



i5m = i^J■L - t^R)/fl, 



(27) 



Which-way information and quantum erasers 



and 



(28) 



Equation ([26| evidentiates that the interference contri- 
bution is an oscillating signal. The phase of the oscil- 
lation, however, depends on time and oscillates with a 
frequency determined by the difference 5^ between the 
chemical potentials, while the phase offset is determined 
by the relative phase between the two quasi-condensate. 
Since the two quantum gases are independent, this rela- 
tive phase will be defined only for a single experimental 
run lliiiE. 



Some remarks on the oscillating behaviour of the in- 
terference term are now in order. In first place, for the 
example considered in Eq. (|26l) the photon flux oscillates 
in time with frequency 5pL. This oscillation is observed, 
however, only when one chooses a time- window at the de- 
tector Ai such that <C 1. Indeed, the photon flux is 
composed by two components: one is constituted by the 
atoms extracted from the left condensate, which interfere 
with the atoms extracted from the right trap during the 
time interval At (flrst term on the RHS of Eq. (gS])). The 
second component is constituted by the atoms extracted 
from the right trap, which interfere with the atoms ex- 
tracted from the left condensate during the time interval 
M (second term on the RHS of Eq. ^^). This interfer- 
ence is analogous to the interference between two lasers 
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at different frequencies: interference fringes are observed 
provided that the time window of the detector is suffi- 
ciently small so not to resolve the detuning between the 
lasers [i^ 11^ • As in the case of two independent lasers, 
there is no well defined phase offset: the relative phase is 
defined only for a single experimental run, while the av- 
erage over a statistically significant number of runs gives 
no interference pattern [S, 43, 44, 46]. 

Differing from the situation of two laser beams p8| . 
however, photon scattering here creates correlations 
among the scattered photon, the corresponding outcou- 
pled atom, and the quantum gas. This correlation is a 
form of a which-way information, that in general washes 
out any interference in the photon flux [50] and can be 
considered as a form of photon-atoms entanglement [sij . 
Oscillations in the photon flux, and hence interference, 
can be recovered if (i) the beam of atoms outcoupled 
from one well overlap spatially with the wave function of 
the atoms in the second trap and (ii) if there are either 
initial correlations between the atoms in the two wells or 
if they both possess a non- vanishing condensate fraction. 

Condition (i) corresponds to a quantum eraser (3^ : 
when it is fulfilled, in fact, the outcoupled atom is no 
more entangled with the scattered photon. We will see 
in the following that it can be fulfilled for certain geome- 
tries and after a certain time, corresponding to the time 
needed for the outcoupled atoms to travel to the second 
well. 

Condition (i) alone is however not sufficient. The pho- 
ton is not only correlated with the outcoupled atom, 
but also with the scattering system. Therefore, only 
when either the atoms in left and right system are corre- 
lated (via e.g. tunneling events which may wash out the 
which-way information [52| ) or when these correlations 
are partly classical, as it occurs when there is a finite 
condensate fraction in both systems, then interference 
can be reestablished. This is essentially condition (ii). 

We remark that oscillations of the photon flux as a 
function of time have been observed in this setup in 
Refs. [s^, Here, they have been discussed in terms 
of Josephson Junction, taking into account the possibil- 
ity that the outcoupled atoms from one condensate can 
be then retrapped in the second well by means of an in- 
verse Raman process. The discussion we just reported 
provides a quantum optical interpretation of the phe- 
nomenon, that is valid for short times, when the quan- 
tum state of the scattering atoms is not substantially 
perturbed by the scattering process. 



2. Onset of the time-dependent oscillations 

We now analyse the onset of oscillations in time. As 
we mentioned, the photon flux starts oscillating after a 
finite time has elapsed from the beginning of the ex- 
periment. [sO, [3l|. This can be also seen in our the- 
ory, by performing the integral over k in Eq. (1261) . 
In App. |X] we show that Fj can be recast in the form 



Fi{t) - FL^Rit) + FR_,L{t), with 



FL^R,{t) w TRee'f'^^*-'^^^) / dre'("<'-^)^/iL7?(d, r) (29a) 



FR^Ut) « rRee-'^''''*-^^^) J dTe'('^<'-"+^^)^/iflL(-d, r) . 



(29b) 

Here, 



h,i{d, t)= j dr/, (r + d - v^t) fi{v) 



(30) 



is the overlap between the left and the right condensate, 
with one being displaced by the amount VqT — d, with 
the recoil velocity 



(31) 



that is acquired by the outcoupled atom by scattering 
the photon. This overlap vanishes at time t = 0, i.e., 
when the outcoupling lasers are switched on (recall that 
the two clouds initially do not overlap). We note that the 
overlap integral is zero at all times if d and q are orthog- 
onal, while it may become maximum after a certain time, 
when d and q are parallel, say, pointing along the posi- 
tive X axis. In this case, the component FL^R{t) may not 
vanish and can be interpreted as the contribution to the 
interference flux Fi[t) from the outcoupled atoms that 
propagate from left to right (for the term FR-^L{t) this 
is just opposite). Interference will set in, provided that 
a sufflciently long time has elapsed to warrant overlap. 
This corresponds to times t > tc with 



tr 



Vn 



where dcs = d — (^l + Cfl)/2 is the effective distance of 
the two systems taking into account their width ^j. In 
terms of an interference experiment, for times t > tc the 
which-way information has been erased and oscillations 
in the photon flux can be observed. 

In interferometric setups the amount of visibility and 
which-way information are related by an inequality 5Q] . 
It is therefore useful to determine a visibility of the oscil- 
lating signal, as it contains information on the properties 
of the scattering systems. A visibility can be defined 
for sufflciently long times, averaging over several oscil- 
lating periods after the instant tc^ and will be propor- 
tional to the amplitude of the oscillations of the photon 
flux, hence to the product of the condensate fractions of 
both systems. We remark, once again, that oscillations 
can be observed only in a single experimental run, while 
they will disappear after performing an ensemble aver- 
age. Therefore, this behaviour can only be measured in 
systems, whose properties are not deeply modified by the 
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outcoupling of atoms. When this is not verified, the pres- 
ence of a condensed fraction in the right system can be 
revealed by performing an ensemble average over a suf- 
ficient large number of experiments, in which the signal 
is monitored for sufficiently short times, warranting that 
the properties of the quantum gas have not been signif- 
icantly modified, and taking the statistical distribution 
of the intensity of the photon fiux at a given instant of 
time. An amplitude can be extracted from the width of 
the distribution by taking into account the finite width 
of the diffraction function. 



IV. HOMODYNE DETECTION OF A 
QUANTUM GAS 

The theory presented so far will be now applied to some 
specific examples. The main idea is to use the setup in 
Fig. [l]in order to determine the mean value of the field 
operator of a quantum gas, using a Bose-Einstein con- 
densate at known temperature as reference system. Such 
setup is a matter-wave analogon of homodyne detection 
in quantum optics. The individual elements can be so 
identified: the BEC acts as a local oscillator, the outcou- 
pling procedure as beam splitter, the relative phase can 
be varied by changing the interwell distance. The infor- 
mation on the atomic gas is carried by both scattered 
photons and outcoupled atoms: homodyne detection of 
the scattered field, hence, allows one in principle to de- 
termine the mean value of the quadrature of the atomic 
field operator. 

A. Interference between two Bose-Einstein 
Condensates 

We first discuss the case in which two BEC are trapped 
in the left and right well, respectively, and are both illu- 
minated by the pump and probe beams. The outcoupled 
atoms from both condensates propagate along the direc- 
tion determined by the distance between the minima of 
both wells, and overlap after the time has elapsed. 
The scattered photons are revealed at a detector in the 
far-field. This setup has been realized in the experiment 
reported in Refs. [sol [3l|. where time-dependent oscilla- 



tions in the atom and photon flux were measured. Here, 
we apply the theoretical model developed so far and find 
it reproduces qualitatively the results of these experi- 
ments. Moreover, we discuss the results in the light of 
question addressed in the present work. 

We assume two BEC with equal number of atoms Nc 
and temperature T , that are confined either in the left 
or right well. The wells are described by the potential 

W = - i-j ) + ^i.L^"^ (32) 

with V{v) — \m{ijj1x'^ + w^t/^ + ^Iz"^), and IS.V denotes 
the constant offset between the two traps. For simplicity 
we take that both BEC are at zero temperature, T = 0, 
and the atoms weakly interact, such that the contribution 
of the noncondensed atoms to the photon flux is small 
and can be neglected. In this limit the chemical potential 
of both condensates is equal and given by /i(0). With 
definition ([5^ then fin = /i(0) and /i^ = m(0) + 

The photon flux is evaluated using the Thomas-Fermi 
approximation for the condensate wave functions [29l |. 

/,(r) = [(M0)-y(r-r,))/5]'/', (33) 

where /Lt(0) = {l5Ncas/a)^^^ huj/2 is the chemical po- 
tential at zero temperature and a = \J hj {mCo) is the 
size of the ground state of a harmonic oscillator with fre- 
quency uj — [wxijJyijJz)^^^ ■ The condensate macroscopic 
wave function has size ^ = yj2^,{Q) / {mu'j). Its 

Fourier transform reads /L(k) = e''' '*/^/o(k) (for the 
right condensate //{(k) — e^''' '*/^/o(k)), where /o(k) is 
the Fourier transform of the macroscopic wave function 
centered at the origin and is real valued. In particular, 

/o(k) = .o^^, (34) 



where kq = y Ibn^Ncrx ry 'ri ' /2 is a scalar and J2(po) 
the Bessel function of second order of the variable po 
deflned as = klr^°^^ + fc^r^"^^ + fc2^i°^^ see 0. The 
integral for the interference contribution to the photon 
fiux in Eq. can be then cast in the form 



Fi{t) = 2ttkIT 



dk |J2(Po)P 
(27r)3 



■Re 



(35) 



This equation can be simplified, taking that both d and q point along the positive x-direction. Following the derivation 
reported in App. [B] we find that the total fiux can be approximated by the expression 



F{t) « 2TTrNcK„ 



1 + Vo cos 6fj,t - ifiLR + (wq - ^) — 8 i - 



d — 2rx 



(36) 



r 



where 



and K{x) = ■\/l/(27r(T2) exp[— .t^/ (2(t^)] is a Gaussian of 



K„, = (i^(wq + A'(wq - r2 + 5^l))/2 
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width cr^ = 2.5(wq/r:c)^. Equation (p6l) shows that the 
photon flux starts oscillating for times t > tc, with tc = 
(d — 2rx)/vq, namely, when the outcoupled atoms from 
one condensate have reached the second one. The time- 
dependent oscillations have frequency 6fj, and amplitude 



Vn = 2 



K{ujq - n) 



K{ujq + if (wq -n + ^^J.) 



(37) 



The expression in Eq. ((36|) has been obtained neglecting 
the momentum dependence of the trapped clouds (this 
approximation is verified in App.[B]) and the contribution 
of the noncondensed fraction. We note that, although os- 
cillations as a function of time are observed provided that 
(5/i 7^ 0, nevertheless their visibility (which corresponds 
to Vo) is always smaller than unity, Vo < 1 for J^u 7^ 0. 
This can be understood considering that due to energy 
conservation in the scattering process, a scattered pho- 
ton carries partial information about from which cloud it 
was scattered. 

We now turn our attention to the phase offset, charac- 
terising the oscillations of the photon flux as a function 
of time in Eq. ([55]) . This phase offset is here given by 
the quantity 6 — ipLR, where ipLB. is the relative phase 
between the two BEC, which can take any value (we refer 
to the discussion in Sec. IIIIC|) . and by the quantity 




FIG. 2: Photon flux F{t) (in units of the background contri- 
bution FB{t)) as a function of time (in units of &fi/2n) and 
of the angle a at which the atoms are emitted. The photon 
flux is computed by numerically integrating Eq. (|35p . The 
condensates are composed by A'^ = 10® sodium atoms in a 
spherical harmomic trap of frequency u} — 325Hz. The scat- 
tering length is as = 55ao with ao Bohr radius. The other 



parameters are — 6^, S)j,o 
a = and (fLR ~ 0. 



27rlO''Hz, d = 5r^, Q 



Q = —{ujq-n + S^) , 



(38) 



which is the phase the outcoupled atoms accumulate 
when travelling from left to the right well (The corre- 
sponding phase offset for the outcoupled atoms travel- 
ling from the right to the left well is given in App. IB|) . 
This phase can be varied by either tuning the laser fre- 
quency, and thus ft, or changing the distance between 
the minima of the two wells. The phase offset in Eq. 
((38)) agrees with the one reported in Ref. [3l|, that was 
derived from a phenomenological model. In Ref. [sij , in 
particular, the rate of change of the difference in phase 
offset between atoms outcoupled to the left and atoms 
outcoupled to the right was measured as a function of 
the Bragg-frequency il, and found in agreement with the 
prediction of Eq. (|B11[) . 

We also note that our theory predicts correctly the 
oscillation period observed in the experiment, and also 
accounts for the fact why the interference current needs 
some time to build up. The linear response treatment 
we apply, however, does not predict any decrease in the 
visibility of the interference pattern with time, in marked 
contrast with the experimental results. Possible reasons 
for the decrease in visibility in the experimental data are 
depletion of the condensate and heating of the sample, 
which could be due to spontaneous Rayleigh scattering 
events [s^l- These effects are not taken into account in 
our model, but could be introduced by means of quantum 
Langevin equations, using a formalism similar to the one 
developed in Ref. 



atoms, thus for geometries where the direction of emis- 
sion q of the outcoupled atoms forms an angle a with 
the vector d. For simplicity we assume that both vectors 
lie in the x — y plane. The calculation is performed by 
numerically integrating Eq. (|35|). 



Figure [2] displays the photon flux F{t) as a function 
of time and of the emission angle a. Oscillations as 
a function of time are observed for values of the angle 
about a — 0, tt, corresponding to the atoms propagat- 
ing in the direction parallel to ±d. The period of the 
oscillation is 27r/(5/i, independent of the angle of emis- 
sion. The oscillations disappear for angles in the interval 
■""/S ^ o. < 7tt/8: the photon flux is here solely given by 
the background contribution FB{t). For these angles, in 
fact, at all times there is no spatial overlap between the 
wavefunction of the outcoupled atoms and the wavefunc- 
tion of the trapped atoms in the other well. Indeed, from 
a simple geometric argument one finds that the overlap 



vanishes for angle a > arctan ■ 



1 and a < TT - 



We finally discuss the dependence of the photon fiux 
on time and on the angle of emission of the outcoupled 



This implies that the setup must be so constructed, that 
the atoms outcoupled from one well could in principle be 
transferred, by a Raman process, into the second well. 
This property is the key element on which the analogy 
to a Josephson Junction has been drawn [sij. It is also 
basically the way in which a quantum eraser is realised in 
this setup. We refer the reader to the discussion on the 
properties of this quantum interference process in Sec. 

Imtl 
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B. Thermometry of a Bose-Einstein condensate 

We now show how the setup of Ref. [s^ could be 
used in order to determine the temperature of a Bose- 
Einstein condensate, using as reference a second BEG 
at known temperature. For this purpose, we make the 
same assumptions as in the previous section, with the 
difference that while the left gas is a BEG at temperature 
T — 0, the gas in the right well is a BEG at temperature T 
to be determined. We further assume that the Thomas- 



Fermi approximation can be performed also for the right 
condensate. Hence, the chemical potential of the second 
condensate can be written as n(T) = n{0){l~T^ /T^)^^^ , 
with Tc the critical temperature for the noninteracting 
gas in a harmonic trap, while the size of the macroscopic 

(T) 

wave function of the right condensate scales with — 
rf\l - T3/r3)i/5^ see [29|. Using these relations the 
integral for the interference contribution to the photon 
flux can be cast in the form 



/d\f 
(2^)3 



dk \J2{pq)J2{pt)\ 



Pt 



xRe 



,i(5M(T)t-k.d-v.i,fl) Jt(^ - CJk+q) + c-H'5M(T)t-k.d-^^«) Jt(^ _ ^^(j,) _ ^^^^^ 



(39) 



withpT — (1 — T^/Tc )^/'^po- The integral depends on the 
temperature both through a scaling factor as well as the 
function J2{pt)/Pt in the integrand, while the phases 
depend on the temperature of the second BEG via the 
chemical potential of the right BEG, which enters in the 
quantity 

dfi{T) = (5^(0) - T{N, T) /h (40) 

with T{N,T) = fi{T) - fi{0). We note that Fi{t) os- 
cillates as a function of time with both frequency and 
amplitude which depend on the temperature of the sec- 
ond condensate. For T <^ Tc, in particular, one finds 
that Fi{t) (X y/nciT), where nc(T) = Nc{T)/N = 
1 — (T/Tc)'^ is the condensate fraction in the right trap. 

Figure ([3^) displays Fj (t) as function of time and tem- 
perature T of the right condensate, when the trap is 
spherical. The oscillations are visible for times t > tc- 
The oscillation frequency depends on T, as one can 
clearly observe from the figure [ssl ]. The dependence of 
the amplitude of the oscillation on the temperature is 
visible in Fig. ([3^) and is singled out in Fig. ([SId), where 
the amplitude 

C(T) = [maxiFiit)) ~ min(F/(i))] , (41) 

evaluated at times t > tc, is displayed as a function of 
T in units of C(0). The red dashed curve represents 
the squared root of the condensate fraction, y/iyjiT), 
and is reported for comparison. Function C(T) decreases 
monotonically as the temperature increases from T = 0, 
and vanishes at the critical temperature T = Tq (which 
takes into account the effect of the interactions and is 
such that nc(T*) = 0) H. 



C. Monitoring the Mott-Insulator/Superfluid 
transition 

We now assume that the quantum system in the right 
trap is a two-dimensional optical lattice, where the pa- 
rameters can be tuned so that the ultracold atomic gas is 
either in a Mott-insulator or superfluid state, while the 
left system is a BEG at T = 0. Within a mean-field treat- 
ment, in the superfluid phase a nonvanishing superfluid 
order parameter would give rise to a time-oscillating be- 
haviour of the photon flux {EE]. Before we start, we recall 
that the Mott-insulator/superfluid transition in two di- 
mensions has been experimentally characterised in Ref. 
[56| . and the condensate fraction has been measured in 
Ref. [13, while Bragg scattering by this system was re- 
cently reported in Ref. [2]| . 

Here, we argue that the photon flux may allow one to 
detect the superfluid order parameter, and thus monitor 
the quantum state of the system across the phase tran- 
sition. We note, however, that a measurement based on 
outcoupling atoms from the optical lattice will change 
the phase of the atoms, if they have been prepared in the 
Mott-insulator phase. For this specific situation, indeed, 
the approximation we make by treating photon scattering 
as a weak perturbation of the many-body state remains 
valid as long as the scattering is performed over a time- 
scale in which the creation of the local defect does not 
affect the rest of the system. This time scale will be in- 
versely proportional to the sound velocity with which a 
perturbation propagates in the corresponding superfluid 
state. This restriction shows that this measurement pro- 
cedure will not provide a reliable estimate of the super- 
fluid order parameter close to the phase transition point. 

Keeping this in mind, let us now assume that the two 
atomic gases are at T = and tightly confined in the x- 
direction by a harmonic potential, so that the motional 
degrees of freedom in the x-direction are frozen out and 
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T 

(b) 

FIG. 3: (color online) (a) Interference term of the Photon flux, 
Flit), Eq. p9[) . (in units of the background current obtained 
when both condensates are at zero temperature Fb{T = 0)), 
as a function of time (in units of 2n /S/j.) and temperature T 
(in units of critical temperature T^). (b) Amplitude C{T), 
Eq. (I41|l (in units of C(0)) as a function of the temperature 
T (in units of T^). The red dashed line corresponds to the 
squared root of the condensate fraction in the right BEG, 
y^ic(T). The other parameters are as in Fig. [2] 

the systems are essentially two dimensional. The left po- 
tential is assumed to be a highly-anisotropic harmonic 
trap and right potential is a two-dimensional optical lat- 
tice. They read 

VUr) ^ ^mcvlx'' + ^mcvliy^ + z^) , (42) 

Vn{r) = Imcjlx' + Vo (sin' ^ + sin' ^] . (43) 
2 V "0 do J 

with uj± ^ LUx, while Vq is the height of the optical lattice 
potential and do the spatial periodicity in the transverse 
directions. 

Assuming that the left quantum gas forms a weakly- 
interacting Bose-Einstein condensate at T = 0, and that 
fi^x ^ Mi(0)7 with /iL(0) the chemical potential of 




r 

(b) 

FIG. 4: (color online) (a) Interference contribution Fi{t) as 
a function of time (in units of 2-k/5^) and of the parameter 
r = zJ/U. Fi{t) is given in units of TM, where M = 50 x 50 
is the number of lattice sites, (b) Amplitude C(r) (in units 
of C(0)) as a function of r. The parameters are Vo — 10&j_r, 
xo = O.lSdo, 5^ = 0.2tJij, d = 20do, |q| = 2-k /do, a;q = ur, 
with ujR = hiv'^ / {2mdo) , and (pLR = 0. Moreover, ^r/U = 
^ - 1 and riB = (&)%=o.2/dg ~ 0.2/dg. 



the left condensate, then the macroscopic wave function 
reads $L(r) = /o(2;)v'«Le"^^ , with rt^ is the planar den- 
sity of condensed atoms in the left trap (which is here 
assumed to be constant, under the condition that the 
laser beams illuminates the center of the condensate) and 
fo{x) = . ^ /(2a;o) is the motional ground state 

Y xo\/ir 

of the harmonic oscillator in the x-direction. The wave 
function for the right quantum gas is here given assum- 
ing the tight-binding limit. Using the Wannier decompo- 
sition, in the single-band limit one finds 

■ipR{r) = fo{x) ^ wi{y)wm{z)hi^m , (44) 
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where wi ) is the Wannier function centered at ^; = Ida 
— y, z) and bi „i annihilates an atom at the lattice site 
with y — Mq and z = mdo- The corresponding Bose- 
Hubbard model reads (SSj 

l.rn 

+ U^ni^m{,ni^rn - ^) - ^J'R^nl^rn (45) 

with n/ „i = bj ^^bijn the number of atoms at site {l,m), 
J the hopping rate, ?7 the onsite interaction, and iifj the 
chemical potential for the right system. We note that the 
two-dimensional limit is consistent when hujx S> /i/j. In 
the superfluid phase the expectation value of operator bim 
over the lattice ground state does not vanish and is a con- 
stant, (6i,m) = (b), and the wavefunction for the super- 
fluid fraction reads $fl(r) = fo{x){b) J2i,m^iiy)wm{z). 

The photon flux, and in particular the interference con- 
tribution, can now be evaluated using these quantities. 
For simplicity, in the following we assume that the geom- 
etry of the laser is such that q is parallel to the vector 
d. An analytical expression can be determined using the 
Gaussian ansatz for the Wannier functions [s^ and shows 
that Fi (t) cx (6) and is a convolution of signals oscillating 
with phase ^lr + kxd — 5^t, weighted by the occupation 
of the momentum kx fulfilling energy conservation. 

We evaluate Fi (t) by numerically integrating Eq. . 
The superfluid order parameter (6) is determined nu- 
merically within the mean-field treatment of the Bose- 
Hubbard model [111 . Figure EJa) displays the interfer- 
ence contribution Fi{t) as a function of time and of the 
ratio r — zJ/U between hopping strength and the onsite 
interaction, with z the coordination number. This ratio 
is assumed to be varied by changing the onsite interac- 
tion U, while the lattice depth Vq is kept constant. For 
small values of r = zJ/U, when the atoms in the right 
well are in the Mott-insulator state, no interference sig- 
nal is found. At a critical value {zJ/U)c (unit filling) [601 
the right system undergoes the quantum phase transition 
to the superfluid state and the interference photon flux 
Fi{t) starts to oscillate in time with finite amplitude. We 
study the amplitude of the oscillations as a function of r 
by characterizing the quantity 

C(r) = max(i^7(i)) - min(F/(t)) , (46) 

where the minimum and maximum of Fi{t) are found 
over a time interval large compared to the inverse os- 
cillation frequency of Fj{t), similarly to Eq. (PT|) . This 
quantity is displayed in Fig. |4l^b) as a function of r. For 
r > Vc the amplitude C{r) is nonvanishing and increases 
proportional to the magnitude of the superfluid order pa- 
rameter. 



V. FIRST-ORDER CORRELATION FUNCTION 

Several methods have been proposed in the literature 
for determining the first-order correlation functions, that 
are based on time-of-flight techniques, see for instance 
Refs. d, HJ [HI]. Measurements of the first-order cor- 
relation function have be en p erformed on Bose-Einstein 
Condensates in Refs. [63l. [64|. 

In the following we will show how an extension of our 
previously considered setup may allow one to determine 
the spatial first-order correlation function of a quantum 
gas by the scattered photons. In this case, the lasers shall 
illuminate two spatially separated regions of a quantum 
gas confined in a single well potential. More specifically, 
the spatial dependence of the laser-atom interaction in 
Eq. |3]will be characterized by the Gaussian envelope 

|7,(r)|=7oexp(-(r-rj)VAr2)/(V^Ar)i/2, 

with width Ar ^ |d| between the two regions (We re- 
mark that the excitation could be realized with subwave- 
length resolution [1, [3, [fil]). The corresponding photon 
fiux then reads 

t 

F{t) = r ReJdt'JdrJ dr' 

j,k=L,R Q 

X7, (r)7fe (r')/(r, t- r',V)G^^^ (r, t; r', i') (47) 

where 

G^'Hr,t;r',t') = {^l{r,t)Mr',t')) (48) 

is the first-order correlation function and f{r,t;r',t') 
is defined in Eq. (jlSp . Let the pump-probe excitation 
be a pulse of mean duration t such that uit ^ 1, but 
LOat ^ 1, with LOa the typical frequency characterising 
the excitation spectrum. In this limit we can approxi- 
mate G(i)(r,t;r',t') ~ G^^) (r, 0; r', 0) in Eq. For 
convenience, we denote by G*^^^(r;r') = G^^-' (r, 0; r', 0) 
the spatial correlation function. For Ar sufficiently small, 
so that it can be approximated with a 5-like excitation, 
the photon fiux can be recast in the form 

F{t) ~ if (l + Re{e"i 'iG(i)(d/2;-d/2)}/no(d/2)) 

(49) 

with if is a constant, determined by the details of the 
excitation scheme, and no(d/2) — G^^^ (d/2; d/2) is the 
density at r = ±d/2 (assuming the system has reflec- 
tion symmetry about r = 0). Therefore, the photon flux 
exhibits oscillations with a visibility which is determined 
by the spatial flrst-order correlation function. By varying 
the scattering wave vector q one would thus measure the 
first-order correlation function as a function of d. Re- 
alistic excitation schemes are of course characterized by 
finite spatial resolution Ar, which results in averaging 
the correlation function in Eq. over the finite size of 



13 



the illuminated region and hence to a diminution of the 
contrast (52|. 

This setup could be extended to measure time- 
dependent correlation function by applying a pair of laser 
pulses: Assuming that the photon scattered after the first 
pulse can interfere with the photon scattered after the 
second pulse, then the photon flux at the detector will 
exhibit oscillations whose amplitude is proportional to 
the first-order correlation function (|48l) . A possible real- 
ization could use a mirror placed in front of the quantum 
system, as realized in [66j . 

VI. CONCLUSIONS 

In this work we have discussed a setup which allows 
one, to measure the condensate fraction and the first- 
order correlation function of a quantum gas by means of 
photo-detection. The photons are scattered by the quan- 
tum gas in a pump-probe type of excitation, such that 
the scattered photon is associated with an outcoupled 
atom, with which it is entangled. In addition, photon 
and atom are correlated with the quantum gas. This 
correlation is detected in the photon flux, provided that 
certain conditions are fulfilled, which we have identified 
and discussed. 

Our analysis is based on the impulse approximation 
[43I , hence on neglecting the back action of the scattering 
process on the quantum gas. It is therefore valid for short 
time transients. From the point of view of collecting a 
sufficient statistics, hence, time-of-flight techniques are a 
more convenient tool. Nevertheless, one could consider to 
modify existing techniques ,^uch as the one successfully 
demonstrated in Refs. [1,13, order to access to the 
same kind of information that the setup here discussed 
provides. 

An interesting outlook is to identify a setup, along the 
lines of the proposal Ref. which permits one to 

perform a quantum-non-demolition measurement of any 
correlations function of of the external degrees of free- 
dom of atomic gases, and more in general, which realizes 



quantum-state transfer between matter and light. This 
would open several interesting perspectives for quantum 
communications [csj . 
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Appendix A 

We recast the interference term, Eq. (jl3bp . in the form 
Fi = FL^iiit) + FR^L{t), with 

t 

= FRee'*^* f dt' f -^e~"^''~^^^'^'~''^ 
J J (27r)'^ 


X J drdr'e'(i-'') ('"-''')$L(r')*$i?(r),(Ala) 

t 

FR^L{t) - FRee-'*^*y" dt' J ^^Q-m-'^u-Sf^Kt-t') 


X J dr dr'e'(i-''> $i?,(r')*$L(r). 

(Alb) 

We now show that the term F]^^R{t) is the 

contribution due to the atoms which are outcoupled and 
propagate from the left to the right (right to the left). 
For this purpose we perform the k integral in Eq. (|A1[) 
and obtain 



fL^J^(^) =FRee'(''''*-'^^«) / dr 



dr / dr'e'('i ('-''')-"^)exp 



i m(r' — r)^ 



fL{v')fR{v), (A2) 



where ^plr is the relative phase of the macroscopic wavefunctions defined in Eq. 
f = r — r', R = (r-l- r')/2, we can rewrite Eq. (jA2[) as 



h 2t 

With the change of variables 



FL^R{t) = FRee'(*'''*-^^«' J dr (^^^) ' e'^""-"'" 



dr / drexp 



1 m 
h2^ 



hq 
r T 

m 



/i(R-|)/;^XR+|)(A3) 



The exponential in the integral over r oscillates very the main contribution to the integral over r comes from 
fast with respect to the wave functions fj{r). Therefore 



14 



Fq = ^T", where the term in the exponential vanishes. 
By means of the saddle-point approximation we take the 
wave functions at the point fg out of the integral. Per- 
forming the f integration using the Fresnel integral (28j 



oo 

J " VM 



(A4) 



we obtain 



z 



X j dr/i (r + d - v^t) fniv) , (A5) 

where Vq is the recoil velocity, Eq. ([5T|) . The calculation 
that leads to Eq. (jA5l) is exact in the limit of homo- 
geneous atomic systems. In such case the momentum 
distribution of the condensate fraction is a Dirac-delta 
function at zero momentum and all outcoupled atoms 
have exactly the same momentum hq. The condensates 
we consider are confined by an external potential and 
have a finite extension in space which leads to a certain 
width dp in their momentum distribution, and thus to a 
spread in the momentum of the outcoupled atoms around 
the mean value hq. By taking the value of the atomic 
wave functions at the point of the stationary phase of 
the exponential in Eq. (jA3l) one neglects this momentum 
width: The saddle-point approximation is thus applica- 
ble if ?i|q| » 6p. Equation (|A5p agrees with the corre- 
sponding expression in Eq. (j29b ,). For completeness we 
also give the contributions of the macroscopic wavefunc- 
tions to the background contribution, Eq. (I13ap . making 
the same approximations as in Eqs. (|A5|) : 



z 

FlU) « FRe j dTe'('^--")^y dr/L(r)/L(r - z;qr) , 



t 

Fnit) « FRe J dTe"^'^-~''+'''^^ J dr /«(r)/fl(r - ^^r) . 



(A6) 

Appendix B 



We derive here approximate expressions of the Raman 
scattering rate for the experimental setup of 's^, which 
lead to Eq. (1551) . Using Eq. ([55)) for the condensate wave- 
functions in Eq. (jA6[) we find 



FUt) 



giVq 



2<(0) 



dz cos 



G{z), 
(Bl) 



with 



z<(x) = min(2, — — —) . 



(B2) 



The length R = {rxryr^Y^^ determines the typical size 
of the condensates and can be written as Gftl 



R=15 



1/5 



(B3) 



The function G{z) in Eq. (|Bip reads 



4/15 




FIG. 5: (color online) Comparison of the integral G{z) as 
given in Eq. (|B4|) (black solid line) with the Gaussian fit 
Eq. (jBSP (red dashed line). 



G{z) = 2 I dx I rdrfl-x^-r^^^^^ 



z/2 



(1 - (x - z)2 - r2) 



(B4) 



A numerical evaluation of G{z) is shown in Fig. ([5]), and 
is here compared with a Gaussian of the form 



Giz) . -e 



-1.25Z-' 



(B5) 



showing that this function provides a good approxima- 
tion of Eq. ((B4)) . Using Eq. (jBS)) in Eq. ((BT|) one gets 

^<(0) 



FL(t) « FA^c^o J dzcos[{iUq-n)zto]e' 





(B6) 



with to = Tx/vq. For times t > tc, with tc (d— 2ra;)/wq, 
such that 2;<(0) = z^{d) = 2 we find from Eq. (IB6p 



FL{t) « rtoNc J dzcos 



(Wq - ft) — 



-1.252-= 



V Stt 



(B7) 
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When performing the integral in Eq. (jB7|) we neglected 
the imaginary part in the Error-functions which one gets 
from the exact integration. This is a good approxima- 
tion for \{ujq — f2)<o| ^ 5, as can be checked by numerical 
evaluation. For values | (wq — fl)to \ > 5 the exponential 
in Eq. (|B7p is negligible compared to the resonant case 
Wq = ri: The oscillating tails of the Error-functions given 
by their imaginary parts will only play a role for parame- 
ters where the outcoupling efficiency vanishes and which 
are thus not relevant to our treatment. Using the same 
argumentation for the other contributions to the photon 
flux we find 



with 







K{x) = W-O-e-T' 

V OTT 

Fe{q,t) = Q{q)e(^t 



d-2r. 



(B9) 
(BIO) 



From Eqs. (IBSp one then obtains Eq. ([55]) for (7 > 0. 
The phase difference between FL^i^{t) and Fj^^L{t) as 
measured in 31] is obtained from Eqs. (jBSb .d) and reads 



^TNcK{ujq - Q.) , 
2nTNcK{u;q-n)Feiq,t) 



X cos 



5lJ.t - ipLR + (Wq - ri) — 



(B8a) 
(B8b) 



c) 



FR^Lit) ~ 2TrTNcK{LUq - n + Sfi)Fe{-q,t) 



dfit - ipLR - (Wq - 17 + dfi)- 



(B8d) 



e = — (2f7-2Wq-(5/l). 



(Bll) 



Figure ^ displays the total photon flux computed 
from Eq. (pS)) (dashed line) and the one calculated from 
Eq. p5|) (solid line) for the experimental parameters of 
[SOj . Quantitative agreement between the two solutions 
is found, showing that neglecting the initial momentum 
distribution is well justified for the considered parame- 
ters. 
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